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ABSTRACT 

The glass transition is described in terms of thermally activated local 
structural rearrangements, the secondary relaxations of the glass phase. The 
interaction between these secondary relaxations leads to a much faster and 
much more dramatic breakdown of the shear modulus than without interaction, 
thus creating the impression of a separate primary process which in reality does 
not exist. The model gives a new view on the fragility and the stretching, two 
puzzling features of the glass transition. 

The mode coupling theory of the glass transition (Gotze and Sjogren 1992) 
postulates a crossover in the flow mechanism at the critical temperature T c to thermally 
activated hopping on the low-temperature side. This postulate has received strong 
support by recent simulations, which showed that the system passes from the saddle 
points of the energy landscape to the minima at exactly this critical temperature (for a 
review see Buchenau 2003). Thus one indeed expects thermally activated flow between 
T c and the calorimetric glass transition temperature T g , where the system freezes. 

However, if one looks at the temperature dependence of the flow process below T c , 
it is dramatically faster than that of a thermally activated process, particularly for the 
so-called fragile glass formers (Bohmer et al 1993). This fragility seems to be linked to a 
strong stretching of the shear stress relaxation, extending over several decades in time. 
Both phenomena have not yet found a generally accepted explanation. 

The present paper proposes to describe the glass transition in terms of the energy 
landscape concept (Goldstein 1968; Johari and Goldstein 1970, 1971; Stillinger 1995), 
taking the interaction between different thermally activated jumps into account. There 
is increasing evidence (Richert 2002) for a heterogeneous energy landscape dynamics. 
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In the glass phase, the energy landscape idea has been successfully exploited in 
three phenomenological models, the tunneling model (Phillips 1981) for the glass 
anomalies below 1 K, the soft-potential model (Parshin 1994) for the crossover to higher 
temperatures and, finally, the ADWP (Asymmetric-Double- Well-Potential, Pollak and 
Pike 1972) or Gilroy-Phillips model (Gilroy and Phillips 1981) for the description of 
the classical relaxation over a wide range of barrier heights, up to barrier heights which 
become relevant for the flow process, i.e. barriers of the order of 35 ksTg. 

As is well known, two relaxation centers at different positions interact via their 
elastic dipole moments. Here, it is proposed to treat this interaction in terms of the 
weakening of the shear modulus by the relaxation, a mean-field approach. 

We describe the dynamical shear response of the glass or undercooled liquid in 
terms of the apparent barrier density function f(V) defined by 

where SG is the reduction of the shear modulus from all relaxations with barrier heights 
between V and V + 5V and G is the infinite frequency shear modulus. The relaxation 
time ry of these relaxation centers is given by the Arrhenius relation 

t v = T exp{V/k B T) (2) 

with r ~ 10 _13 s. 

One next needs the relation between the true barrier density fo(V) and the apparent 
density f(V). Consider a constant shear strain e applied at time 0. The stress a after 
the time t is given by 

a(t) = G(t)eo = Ge / exp(-t/r v )f(V)dV. (3) 
Jo 

Here Vm is the Maxwell barrier corresponding to the terminal relaxation time tm of the 
shear modulus, the Maxwell time. 

A given relaxation center will tend to jump at its relaxation time ry. At the time 
Ty, the stress is diminished by the factor Gv/G, where Gv = G{jv)- The local distortion 
at the center is assumed to be unchanged. This implies an unchanged contribution of 
the center to the free energy. Thus the reduction of the stress energy by the center is 
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unchanged. But we reduce a smaller stress energy, so f(V) must be larger than fo(V) by 
the square of the factor G/Gy. In a physical picture, the enhancement is due to induced 
jumps of lower-barrier entities in the neighbourhood, which reduce the resistance of the 
viscoelastic medium surrounding the center. 
Gy is given by 

G V = G exp(-r v /r v )f(v)dv » G / f(v)dv, (4) 
Jo Jv 

because the double-exponential cutoff is practically a step function. Therefore it is 



reasonable to define the true barrier density function by 

n 2 



fo(V) = f(V) 



V M 

f(v)dv 



(5) 



since the integral on the right hand side is close to GyjG. 
Solving this integral equation for f(V) one finds 

f(V) = . V M *\ W r (6) 
3 J?" /„(«)*] 

This solution has two important properties. The first is the 1/3- rule: 

/ fa(V)dV=~. (7) 
Jo 3 

The rigidity vanishes when the integral over the noninteracting relaxation entities is 1/3, 
i.e. when the noninteracting secondary processes would have reduced the initial shear 
modulus at time zero to 2/3 of its value. The 1/3-rule has been derived independently 
from energy considerations (Buchenau 2003). 

The second important property is that the breakdown occurs in a rather dramatic 
way, because the relaxing entities at the critical Maxwell barrier value receive a strong 
enhancement, to such an extent that one is tempted to assume a separate a-process 
which has nothing to do with the secondary relaxations. In fact, this more or less 
unconscious assumption underlies most of the present attempts to understand the 
glass transition (Ediger, Angell and Nagel 1996). The above treatment shows such an 
assumption to be unnecessary; what one sees at the glass transition are simple Arrhenius 
relaxations of no particularly large number density, blown up to impressive size by the 
small denominator of eq. (jHJ). 
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One can calculate the barrier density function from dynamical mechanical 

measurements in the glass phase using the approximation (Buchenau 2001) 

2 G" 



f(y,T) 



(8) 



7r GksT 

where V is calculated from the condition ury = 1. 

Though the potential parameters freeze at T g , fo(V,T) will still depend on 
temperature, even in the glass phase. This temperature dependence can be derived 
from a consideration on the free energy of an asymmetric double-well potential (Fig. 
1). In the simplest possible approximation, the free energy of the relaxing entity is 
(Buchenau 2001) 

A 



-k B T\n 



2 cosh 



2k B T 



(9) 



The Boltzmann factor exp(— F/ksT) at T g supplies the asymmetry distribution 

P(A) - cosh(A/2A; B T 9 ), (10) 

so the number density of relaxing entities does indeed increase with increasing A; it 
is not really constant. If one integrates over the asymmetries at a lower temperature, 
one finds a diminution of the effective f (V,T) by a factor of up to 7r/2 in the low 
temperature limit. A useful approximation below T g is 



f(V,T 9 



2arctan (e 



V/2k B T^ 



-{T g /Tf\n{n/2) 



11) 

f(V,T) [ 2 J " ) 

This factor has to be taken into account to calculate f(V, T g ) from measurements of the 

mechanical loss in the glass phase. Once /(V) is known, /o(V) can be obtained from 

eq. ©. 




Figure 1. Asymmetric double- well potential. 
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Fig. 2 (a) shows f (V,T g ) for polystyrene determined in this way from six torsion 
pendulum data at 0.5, 0.64, 0.9, 1, 6 and 10 Hz (Illers and Jenckel 1958; Schwarzl 1990a; 
Schmieder and Wolf 1953; Schwarzl 1990b; Sinnott 1959; Hartwig and Schwarz 1968). 
The data were fitted by the six-parameter expression 

i^) = w^^ v,v " +a ^" ltv ^- (12) 

Two of the parameters, the energy W and the small number C, can be taken from soft- 
potential fits of the low-temperature anomalies (Ramos and Buchenau 1997), because 
the first term is the barrier density of the soft-potential model (Parshin 1994), with an 
exponential cutoff of the soft-potential expectation at a limiting barrier Vq. The second 
term is just an appropriate gaussian with amplitude a± and width parameter gi centered 
at the position p 1 . The results of the fit are compiled in Table I. 
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Figure 2. (a) fo(V, T g ) for polystyrene from torsion pendulum data, together with 
the corresponding f(V,T g ) (b) measured and calculated dynamical shear data at the 
glass transition (references see text). 



Table I: Fit parameters for f (V,T g ). 
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substance 


C 


W 


V 




9i 


Pi 




10" 4 


meV 


eV 


eV- 1 


eV~ 2 


eV 


a — Si0 2 


2.65 


0.34 


0.04 


0.76 


0.54 


5.45 


polystyrene 


7.1 


0.159 


0.1 


0.19 


4.0 


1.0 



With these parameters, the integral over f (V,T g ) extrapolates to 1/3 at Vm = 
l.leV, consistent with the experimental T g of 373 K, defining T g as the temperature 
where the terminal relaxation time r t reaches 100 s. This is the first time a glass 
temperature could be calculated from glass data alone (it is true that the glass 
temperature enters already in the calculation, but one can iterate, and the procedure 
converges quickly). 

Even more impressive, the same parameter combination provides a very reasonable 
fit of the measured (Donth et al 1996) G'{oj) and G"{ui) at the glass temperature (for this 
sample at 363 K). This is shown in Fig. 2 (b), calculated with G = 1.55 GPa, the value 
obtained from the Brillouin transverse sound velocity (Strube, private communication) 
of 1071 m/s, taking the weakening by fo(V) at the Brillouin frequency into account. 
The sharp barrier cutoff at Vm of the model seems to be too sharp for this polymer case, 
possibly because of the crossover from enthalpic to entropic restoring forces. 

Fig. 3 (a) shows the barrier density function in the strong glass former S2O2. 
From measurements at 90 and 180 kHz (Cahill and van Cleve 1989; Keil, Kasper and 
Hunklinger 1993), one finds again the soft-potential model expectation f (V,T g ) = 
2C yjy 3 / 4 ^/ 1 / 4 at very low barriers, with C and W taken from fits to the low- 
temperature thermal conductivity and specific heat (Ramos and Buchenau 1997). There 
is a pronounced exponential cutoff over more than four decades leading to very low 
mechanical damping at and above room temperature, the lowest observed in any glass 
so far. Data close to T g are rather scarce. Therefore it turned out to be necessary to 
include the data (Mills 1974) for G'(uj) and G"(oj) at T g into the fit for the high-barrier 
gaussian in eq. (JI2J)- Thus, one cannot predict T g from glass phase data alone as in 
polystyrene. The fulfillment of the 1/3-rule required G = 30.2 GPa for the fit of Fig. 
3 (b), while transverse wave Brillouin data (Bucaro and Dardy 1974) suggest a value 
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of 33.7 GPa. Nevertheless, having fitted fo{V) both at the low- and high-barrier end, 
it is gratifying to see the good agreement of the very low damping at and above room 
temperature for 37 kHz (Marx and Sivertsen 1953) and 1.5 MHz (Fraser 1970) in Fig. 
3 (a), because these data were not used to obtain the fit. Also, Fig. 3 (b) shows that in 
this case the sharp cutoff gives a much better fit of the stretching than in polystyrene. 

Note that not all glass formers can be fitted by the simple form of eq. (|12p: many 
glass formers show several separate secondary relaxation peaks. 

Finally, let us turn to the fit of the fragility. This requires the temperature 
dependence of a whole function, namely fo(V, T) above T g . Here, the estimate for 
this dependence is based on the Boltzmann factor of eq. ([10)1 . Obviously, the rise of the 
Boltzmann factor with increasing asymmetry A cannot go on forever. We assume that 
it stops at A = b a V, where b Q is a fit parameter of order 1. One expects a return to a 




Figure 3. (a) fo(V,T g ) for S1O2 from measurements at four frequencies and from 
the dynamical glass transition data (b) measured and calculated dynamical shear data 
at the glass transition (references see text). 
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kind of low-energy glass ground state when the asymmetry gets as large as the barrier 
height. 

This assumption allows to calculate a partition function for a given barrier height. 
From this partition function, f Q (V,T) should follow the equation 
fo(V,T) T g smh(b V/2k B T g ) 

fo{V,T g ) T smh(b V/2k B T) 1 ' 

for temperatures above T g . 

With this equation, one can calculate the shift factor = tm(T)/tm(T 9 ) from 
the 1/3-rule, eq. (|7|). The result is compared in Fig. 4 to measurements (Mills 1974; 
Bucaro and Dardy 1977; Plazek and O'Rourke 1971). The fit parameter b Q was 1.4 for 
polystyrene and 0.33 for silica. The difference in fragility does not only stem from this 
difference, but also from the large difference in /o(Vm) at T g : a rise in the number of 
relaxing entities shifts the Maxwell barrier much farther in polystyrene than in silica. 

This fragility fit may be too simple-minded, but it points the way to a new 
understanding. If f (V,T) increases with increasing temperature, the Maxwell barrier 
will decrease. Once one accepts the idea of a relatively high free energy of a symmetric 
double-well potential in the undercooled liquid, the fragility ceases to be a puzzling 
property. 




T /T 
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Figure 4. Angell plot for polystyrene and Si02- The dashed line shows the Arrhenius 
law, the continuous lines the model fits explained in the text. 
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